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Abstract 

The statistical properties of the quantum chaotic spectra have been studied, 

so far, only up to the second order correlation effects. The numerical as well 

as the analytical evidence that random matrix theory can successfully model 

the spectral fluctuatations of these systems is available only up to this order. 

For a complete understanding of spectral properties it is highly desirable to 

study the higher order spectral correlations. This will also inform us about the 

limitations of random matrix theory in modelling the properties of quantum 

chaotic systems. Our main purpose in this paper is to carry out this study 

by a semiclassical calculation for the quantum maps; however results are also 

valid for time- independent systems. 
PACS numbers: 05.45+b, 03.65 sq, 05.40+j 
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I. INTRODUCTION 



In generic hamiltonian systems with many degrees of freedom, the classical dynamics 
shows an enormous richness in structure, increasing with the interaction between degrees 
of freedom. The classical motion is mainly of two types, integrable and chaotic; see [1] for 
details. This paper deals with the quantum properties of hamiltonians whose classical limit 
is chaotic. 

The strongly chaotic nature of underlying classical dynamics intutively suggest us to 
expect some kind of random behaviour in quantum dynamics too. This is because the 
classical dynamics is indeed a limit (h = 0) of quantum dynamics and therefore the nature of 
former should be somehow reflected in the latter. In fact, various analytical and numerical 
studies (see [2] and references there in) have confirmed that the manifestation of chaotic 
behaviour in quantum dynamics occurs through randomization (partial or full) of matrices 
of associated quantum operators. The spectral and strength fluctuations of these operators 
can be well-modeled (up to second order correlations) by one of the various universality 
classes of random matrices. Most common among these are the Gaussian orthogonal and 
unitary ensembles, GOE and GUE, the analogous circular ensemble models COE and CUE 
[3]. The former pertain to autonomous systems whereas the latter have application in the 
study of non-autonomous systems such as quantum maps. 

The presence of RMT type spectra in quantum chaotic systems can be explained by the 
Gutzwiller-semiclassical quantization scheme [4] for time-independent systems which uses 
the elegant technique of path integral sum given by Feynman and relates the chaotic mani- 
folds of classical dynamics to the eigenfunctions of quantum dynamics. A similar formulation 
is given for time-evolution operators of quantum maps too [5]. The spectral fluctuation mea- 
sures can then be determined approximately by using the principle of uniformity [6] which is 
based on the uniform distribution of periodic orbits at large time scales and gives a technique 
to evaluate the sum of periodic orbit contributions. Using this technique for autonomous 
hamiltonians, Berry [7] provided an explicit expression for the semiclassical form factor 
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K 2 (t) -the fourier transform of the 2- level spectral correlation function- for values of r in the 
range r << 1 (the time measured in units of 2irhd, where d is the mean spectral density). 
This result has an exact analogy with the corresponding RMT behaviour; following essen- 
tially the same technique as used by Berry for autonomous hamiltonians, this analogy can 
be proved for quantum maps too [2]. In the region r >> 1 too, the limiting behaviour was 
analysed by Berry using a semiclassical sum-rule which makes use of the properties of the 
function related to the quantum-mechanical density of states. 

Notwithstanding the good agreement between RMT and statistical quantum chaos up to 
2 nd order correlations (long and very long), still there is no reason to believe that the Random 
Matrix Theory (RMT) can model all n th order spectral as well as strength correlations. The 
numerical studies for many systems (e.g Baker map [8], Quantum Kicked Rotor [2], etc.) 
have already indicated that even second order correlation effects, when considered on short 
time scales (i.e very long range correlations) , do not follow random matrix prediction and are 
non- universal. This is the range where the classical dynamics is still diffusive and periodic 
orbits are not yet uniformly distributed. The deviation from RMT in this range agrees well 
with our intution as one should expect RMT to be applicable only on those time-scales where 
the variables associated with classical dynamics are random enough to fully randomize the 
matrices associated with corresponding quantum operators. Moreover, the sum rules for 
the matrix elements of quantum chaotic operators [9] have already been found, differing 
from those of RMT. But a study of higher order correlations between zeros of Reimann-zeta 
function shows a good agreement with RMT [10,11] 

Thus it is relevant to know as to what properties and up to what order the behaviour 
of quantum operators can be modeled by RMT and when it ultimately breaks down. Our 
attempt, in this paper, is to make a comparative study of one such property, namely the 
n th order spectral correlation as all n-level spectral fluctuation measures can be expressed in 
its terms. We fufill this goal by carrying out a semiclassical study of the fourier transform 
of n-level correlation function R n ; the reason to consider the fourier tranform lies in the 
covenience of its analytical as well as numerical calculability. We proceed as follows. 
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The Gutzwiller formulation gives us the density of states as a sum over periodic orbits 
and this gives rise to periodic orbit- interaction terms in n-level density correlation function. 
Berry, in order to obtain result for 2-level form factor, neglected the contribution from these 
interacting terms as a first order approximation (the so called diagonal approximation). 
But, for a complete evaluation of form factor, one has to calculate the contribution due 
to interacting terms. The lack of the knowledge of action correlations handicaps us from 
doing so. One attempt in this direction was made in ref. [12], in which by assuming the 
complete validity of random matrix theory, the periodic orbit correlations were calculated 
from the RMT form factor. These when compared with numerically obtained correlations 
(for Baker map, Hyperbola billiard and perturbed Cat map) showed a good agreement. The 
numerical study of these actions also indicated the presence of an uncorrelated component 
(exponentially larger than correlated part); in this paper we use this fact. We assume that, 
on long time scales as a first order approximation, actions are uncorrelated and calculate 
the n th order form factor. Under this approximation, the result turns out to be same as 
that of RMT which is also confirmed from the numerical analysis for at least two higher 
order fluctuation measures, given in this paper. For a complete calculation of n th order form 
factor, the action correlations should also be taken into account which will determine the 
higher order terms. 

In this paper, we present our semiclassical study for quantum maps but the method 
can easily be generalized for time- independent systems and one obtains the same results. 
Similarly the RMT is given only for Circular ensembles but, once again, the final results 
are valid for Gaussian ensembles too which follows due to GE-CE equivalence for large 
dimensions) . 

This paper is organized as follows: 

In section [11 A| , we briefly review the definition of various random matrix ensembles. 
For later use, we also discuss the relation between the n-level form factor and correlation 
functions. The section|IT"B| deals with a brief review of the fundamentals of quantum maps 



and the earlier obtained results for 2-level form factor. Both the sectionllA and sectionllB 



are included in this paper so as to clarify the ideas used in section|ITI| which deals with the 
higher order correlations and form factos for the quantum spectra with exact symmetry. In 
section|TV|, we numerically study the higher order fluctuation measures, namely, skewness 
and excess for a prototype quantum chaotic system, that is, Kicked Rotor and compare 
them with those of RMT. We summarize our results in sectionjV]. 



II. PRELIMINARIES 

A. Random Matrix Results 

We briefly outline here the results for n th order form factor (i.e fourier transform of n- 
point density correlation function) for the eigenvalues of the equilibrium circular ensembles 
of the random matrix theory. 



1. The Circular ensembles 

The circular type equilibrium ensembles are ensembles of unitary matrices Up in which 
the distinct non-zero matrix elements of U are distributed independently as zero-centered 
random variables; (3 defines the number of independent components of the matrix elements 
of U. There are three such ensembles, characterized by j3, namely COE, CUE and CSE for 
(3=1,2 and 4 respectively. These Universality classes are determined by the invariance of 
the system under time-reversal transformation (TR) (or more generally antiunitary trans- 
formation) and are described by the invariance of the ensemble measure: invariance under 
orthogonal or symplectic transformations for TR-invariant systems and under unitary trans- 
formations for TR-non-invariant ones. The invariance restricts the allowed space of matrices, 
for example, to that of symmetric unitary matrices for orthogonal invariance. 
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2. n-point Correlators 



The afore- mentioned unitarity of U implies that its eigenvalues exp(iEj) lie on the unit 
circle in the complex plane, where exponents E^s are termed as eigen-angles. The density 
of states is then defined by 

N oo 

P( E ) = E E KE- 2tt* - Ej) (2.1) 

j=l k= — oo 
1 oo 

= — V exp{inE)Tr{U n ) (2.2) 

27r n=-oo 

and has the mean value < p >= N/2tt. 

For analytical studies of the spectrum, it is a usual practice to calculate the level density 
correlations. For cases where the level density p(E) can be written as the sum of a smooth 
part < p{E) > and a fluctuating component Sp(E), it is preferable to study the correlations 
R k between the fluctuating parts of the density. The R k s can be defined as follows, 

p(p ^ < 5p(E l )5p(E 2 )....5p(E k ) > 

Rk{El > Ek) ~ < P (E 1 >...<p(E k )> (2 ' 3) 

here 5p{E) = p{E)- < p(E) >, E j = E + IjD for j = 1,2.., k - 1 and E k = E with D 
as the mean spacing and <> implying the averaging over variable E for ranges containing 
sufficient number of mean spacings. By substracting < p(E) > from eq.(2.1) and using 
notation Tr(U n ) = t n , Sp(E) can further be written as follows, 

1 oo 

5p(E) = — Y t n exp{inE) (2.4) 

Z7T n=-oo,^0 

The substitution of eq.(2.4) in eq.(2.3) gives 
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k \ k-l 



Rk(E u ..,E k ) = ^J2 <t n t j2-t 3k > (exp[i(2 j m )E] ) exp[iD(J2 

J \ m=l I E rn=l 

here J2j implies the summation over all indices ji,j2,--,jfc with each index varying from 
— oo to oo except zero (that is, none of the indices take value zero). The averaging over E 
reduces the eq.(2.5) in following form, 
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1 / k \ fe_1 

R k (E u .., E k ) = < tA-tjk > 6 ( £ -M e M iD ( Yl 3ml m )\ 

iV J \m=l J m=l 

Y k-l 



(2.6) 

(2.7) 



here X) implies the J2j subjected to condition that Y^ZtiJm 0- 111 semiclassical analysis, 
instead of dealing directly with R k) it is easier to calculate the k th order form factor, defined 
as follows, 



Kk(n, ..,T fc _i) = J 



exp 



exp 



2m^(E j - E k )r. 

k-l 

2m h T i 

3=1 



Rk(Ei, E k )dEx...dE k _x 
Rk(h, •••) l k -i)dh---dlk-i 



(2.8) 
(2.9) 



Substitution of eq.(2.5) in eq.(2.7) gives K k in terms of the traces, 
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k-l 



E'< *ii^-*-(EX=ii™) > n ^ + ^) 



(2.10) 

J =x- • m=1 1\ 

As in this study we confine ourselves to calculation of K k only for \r m \ < 1, m = 1, 2, ..(A; — 
1), for these values of |r m |'s only those terms in the J2j contribute which have indices 
{ji,j2,-;jk} much less than N. Therefore, the J2'j m eq.(2.5) can be replaced by J^'g i n 
which the indices vary from some value — n to n where n < N with all other conditions 
same. 

The above result for K k can further be simplified by using the recently obtained result 
for the statistics of the traces [13] which indicates that the first few traces ti, t 2 , .., t n of large 
unitary matrices (N large) taken from any of the circular ensembles display no noticeble 
correlation. The ensemble average of each of the traces vanishes (that is, < t n >= 0) for 
all the circular ensembles, due to uniformity of the distribution of E^s, the eigenangles; see 
ref.[13] for details. But note that t n and i_ n are not independent of each other and one can 
show that [3, 13] 



< \t n \ 2 >= (3n/N n«N 



(2.11) 



where j3 is 1 or 2 depending on whether the ensemble is COE or CUE respectively. Therefore, 



in general, < YYj=i tj >— if at least one tj is such that its opposite t-\j\ is not present in 
the product. This product exists only if the following condition is satisfied, 

<n*i> = <nN 2 >=n<w 2 > (2.12) 

Now, as can be seen from eq.(2.10), for A;-odd, every product appearing in the sum, contains 
odd number of t/s and therefore above condition can never be satisfied. This gives, for 
fc-odd, 

Kk-oddin, r 2 , r fc _i) ~ 0, \Ti\ i=l ^ k < 1 (2.13) 

on the other hand, the application of the condition (2.10) gives the following result for k 
even, 

K k -even(ri,T 2 , ...,T k -l) ^ « [ S ( T Pi + T P2) S ( T P3 + T P4 ) • • • ^ ( r Pfc-3 + r Pfc- 2 ) 

P 

\ t pi\\ T P3 1 \ T P5 1 •• • \ T Pk-i I (2-14) 

where a — 1 for CUE and 2 for COE. The p implies the sum over all possible permutations 
of indices Pi,P2, ■■■,Pk-i over the set 1, 2, k — 1. 

The result given by eq.(2.13) and (2.14) are valid only when each \t\ < 1. For cases with 
\t\ ~ 1 or > 1, (i.e n ~ N) one has to take into account the correlation between traces. 
Furthermore though the method adopted here for the derivation of K k result is applicable 
only for ensembles of unitary matrices but the final results are valid also for Gaussian 
ensembles. This follows due to the equivalence of fluctuation measures of the circular and 
Gaussian ensemble in large dimensionality limit. 



B. Quantum Map Vs Classical Map 

A classical map can be described by a canonical mapping M of coordiante variable q and 
momenta variable p at discrete time-step t n to those at t n+1 . 




(2.15) 
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with W(q n+ i, q n ) as the generator of the map such that 

dW(q n+1 ,q n ) dW(q n+1 ,q n ) 

Pn = ~ , Pn+i = ~ (2.16) 

dq n oq n+1 

The nature of the time-step considered can give rise to different kind of maps [1]. For 
example, for time-periodic hamiltonians, it is easier to study the dynamics in terms of fixed 
time steps (i.e the period of the hamiltonian), the related mapping known as stroboscopic 
mapping. For time-independent systems sometimes it is sufficient to consider only those 
steps of dynamics which occur on a definite plane that is, intersections of trajectory with a 
plane (instead of equal time-steps), known as Poincare mapping. 

The quantization of a two-dimensional classical map when the phase space upon which it 
acts is compact leads to the construction of unitary matrices U of a finite dimension N, and 
their semiclassical limit is obtained for iV — > oo. For example, for a canonical mapping on a 
two-dimensional torus (here taken to be a two-dimensional phase space with periodicities Q 
and P in q- and p-directions respetively), the corresponding quantum propagator acts in an 
N- dimensional Hilbert space and is represented by an N x N unitary matrix U. This follows 
as the number of states N allowed to be associated, by quantization, with the finite classical 
space is restricted (Uncertainty Principle); N is determined by the following relation, 

2ixTiN = QxP (2.17) 

here N plays the role of the inverse of Planck's constant, with iV — > oo as semiclassical limit. 

C. Semiclassical Form Factor for Quantum Maps : Symmetry Preserving Cases 

For quantum maps acting in a finite Hilbert space, eqs.(2.3) and (2.9) can be used to 
write the quantum- mechanical two- level form factor K 2 {j), 

K ^)=%£ dr ( 6 p{ E+r ^) 6 p{ E - r ^) ex ^ 2mrT ^) E ( 2 - 18 ) 

Using eq.(2.4) in eq.(2.18), K 2 (k) can further be reduced in the following form, 
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i JV-1AT-1 

Mr) = ^EE exp[ik(E n - E m )\ - NS n>0 (2.19) 

n=0 m=0 

= I|Tr(C/")r-^ n , (2.20) 

Where n = Nr. Now the semiclassical expression of K 2 (r) can be obtained from above 
equation by using semiclassical form of Tr(U n ) which can be expressed as a sum over periodic 
orbits in classical phase-space (r = (q,p)) [5], 

^)=ffEE^W[i^-iW2] (2-21) 

here the amplitude A,- (=rij]9(r n — ro)/^^!" 1 ^ =nj(sinha:j)~ 1 ) of the contribution from 
each (multiply traversed) periodic orbit j, with period n, depends on the stability aj of the 
orbit; for long periodic orbits Aj can be approximated as Aj ~ n,,exp(— etj) = nj-exp(— jrij) 
with 7 as entropy of the classical motion. Wj is the action for one traversal of the orbit, rrij 
is the number of traversals, rij = n/rnj is the period of the orbit with single forward traversal 
and Vj is the Maslov index. The index g refers to the number of symmetric analogs, existing 
in the phase-space, for the periodic orbit. While considering the long-range correlations 
which are mainly affected by the long periodic orbits, it is sufficient to consider \m,j\ = 1 
(and therefore rij = n); this follows due to the principle of uniformity which states that on 
large time scales periodic orbits tend to distribute uniformly in phase space, their density 
increasing exponentially while intensity decreasing. Thus, on large time-scales, long periodic 
orbits which are almost all primitive dominate the phase space. 

The evaluation of \Tr(U n )\ 2 , in semiclassical limit iV — > oo, can be done as follows. As 
obvious from eq.(2.21), \Tr(U n )\ 2 contains terms of the type exp[i(Wj — Wk)/K\ and there- 
fore, for a complete evaluation of K 2 (t), it becomes important to study the distribution of 
amplitudes Aj and actions Wj. But, in semiclassical limit H — > 0, the significant contribu- 
tions comes only from those orbit-interactions for which Wj — Wi< o(h). The contributions 
from other orbit-interactions becomes negligible, in this limit, due to the presence of the 
rapid oscillations leading to destructive interferences. Thus, for the leading order semiclassi- 
cal asymptotics of \Tr(U n )\ 2 , one needs to consider only "diagonal" terms [7] with Wi ~ Wj 
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which, in large-n limit (such that n/N = r << 1) can be evaluated by invoking the Hannay's 
sum-rule for the amplitudes [6], 

\Tr{U n )\ 2 ^g 2 Y.^ (2-22) 

3 

Now by using Hannay's sum rule for the intensities, which comes from the principle of 
uniformity [5] and is given by 

E^(M-^)=^ M«l (2.23) 
one can obtain the two level form factor K 2 (t) [7,2] which turns out to be same as that for 
random matrix ensembles, under small- r approximation (eq.(2.14)), 

K 2 (t) ~ g\r\ (2.24) 

Note that above result is valid only for \r\ « 1 i.e \n\ « N. This limit of validity comes 
into existence due to considerations of only diagonal terms in the evaluation of \Tr(U n )\ 2 . 
For cases |t| ~ 1, one needs to consider the contributions due to the constructive interfer- 
ence of very long periodic orbits (with Wi — Wj ~ o(h)) too which once again requires an 
understanding of the distribution of periodic orbit actions. Moreover, in the derivation of 
spectral densty in terms of periodic orbits, the quasi-energy is assumed to be complex with 
a very small imaginary part e (required to avoid the divergence of the formula, occuring for 
real energies). Due to the finiteness of e, the periodic orbits with period n > n* (that is 
the oscillations with energies 5E < e where 5E = Ti/n and e = h/n*) can not be taken into 
account in this formulation. 

The result obtained in eq.(2.12) for semiclassical form factor is same as that for exact 
symmetry classes of RMT, under the same limit. For example, g = 2 and 1 give correspond- 
ing COE and CUE results respectively [2]. 

III. HIGHER ORDER SPECTRAL CORRELATIONS 



In this paper, we restrain ourselves to the study of the /c-level correlation functions in 
short time limits (where period n of the longest periodic orbit, in phase space, is much 
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greater than unity but much less than N i.e 1 << n « N or \t\ = n/N « 1 ) which 
will allow us to use principle of uniformity in the evaluation of cross multiplication of k 
periodic orbit contributions, thus simplifying the calculations. For simplicity and to explain 
our method, we first calculate the 3 rd and A th order correlations and then generalize it to 
k th order correlations. 



A. 3 rd -Order Correlation 



For simplicity, let us first calculate the 3 rd order correlation function. The substitution 
of eq.(2.4) and eq.(2.21) in eq.(2.3), with k = 3, gives us 

R 3 (E u E 2 , E 3 ) = |L (]T ]T AA,A k 

\ ijk mi,mj,m k =±l 

( ( „ 2tt Wi\\ ( ( „ n 2ir Wj 

exp \ irrii [riiE + n^— + — I I exp \ irrij [rijE + n^— + 



exp ^im fe + 



(3.1) 



Here <> implies a local averaging with respect to E, that is, the energy averaging over 
ranges which are classically small but quantum mechanically large so that a large number 
of levels are included. For example, a good choice is to take the size of the averaging range 
to be E itself i.e to define 



<f(E) > E =±J q E f(E>)dE> 



(3.2) 



Therefore the variation of amplitude, in the above equation, with respect to energy is very 
small (amplitude being a classical quantity) and can be ignored. This gives 

R 3 (E 1 ,E 2 ,E 3 ) = — X! AiAjA k exp (77^4 + ^71^2) 

ijk mi,m,j,m k =±l 



N 3 



N 



(exp [(mini + rrijUj + m k n k ) iE}} ( exp (rriiWi + rnjWj + m k W k ) 



(3.3) 



Due to averaging over E, the contribution of various terms, in eq.(3.3), will be determined 
by the fact whether their exponents contain E or not; the terms containing a factor of type 
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exp[iEn] will not make any contribution. Thus we can divide all the terms into following 
two classes. 

Case (1): Terms with allrii,nj,nk of the same sign (i.e either all positive or all negative) 

On averaging over E, the contribution of these terms to R 3 turns out to be zero due to 
presence of a factor of type exp [±i (\rii\ + \rij\ + \rik\) E] 

Case (2) Terms with any two among (rii,nj,nk) with same sign (+ or —) and third one 
with opposite sign. 

The terms under this case contain a factor exp [± (rij + rij — rik) iE] (and its permuta- 
tions). As mentioned above, these terms will make a non zero contribution iff n k = rij + rij 

(or Hi = rij + n k , rij = + n k ). Thus eq.(3.3) can be reduced to the following form 

3 

R 3 (E u E 2 ,E 3 ) =j^Yl A i A j A k 

ijk 

£ exp {nA + n,-€ 2 )) < exp (W t + Wj - W k )) > 

\ ( 1771 \ 

{mi x - nj £ 2 )j < exp ^— (Wi - W 3 + W fc )J > 
+ exp ( nj £ 2 - ndi)) < exp (Wj -W i + W k )^j > (3.4) 

with 2 nd and 3 rd terms corresponding to rii = rij + n k and rij = rii + n k respectively. 

To evaluate terms of type < e^ Wi+Wj ~ Wk ^ h > we proceed as follows. Here W i: the action 
of a periodic orbit with period n«, can also be written as a sum of Hi single step actions 
Wi = Y^Lq 1 Wi(qi+i, qi)\q n = qo - For strongly chaotic dynanics and on large tine scales, these 
single step actions can be regarded as independent variables with a pair-correlation coefficient 
decaying exponentially to zero. An extension of central limit theorem therefore implies that 
Wi's are Gaussian random variable on large time-scales. Hence 9 = Wi + Wj — W k will also 
be a Gaussian random variable with mean zero, the variance of which is given as follows, 

Var6 = < (Wi + Wj - W k f > (3.5) 
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< W? > + < Wf > + < Wl >~ 3T 



(3.6) 



Here T is the average time-period of periodic orbits given by T = h/6E where 5E is the 
energy range over which average is taken. 

As on large time-scales, the phase-space is densely and uniforemly covered by periodic 
orbits and a typical trajectory can be approximated by a very long periodic orbit. This 
permits us to approximate the average of exp(iWj) over all periodic orbits by a phase- 
space average. This gives < exp(z#) >= exp(— VarO) = exp(— 3T). Here, in eq.(3.6), the 
correlations between various actions i.e terms of type < WiWj > has been approximated to 
zero (Wj assumed to be random variable). But, as mentioned in section|, the correlation 
between actions are not entirely zero, that is, W^s are not exactly random variables. The 
distribution function P{6) of these actions can be written as P{6) = P ra ndom + Pcorreiated 
where random part of the distribution P ran dom dominates the non-random part P CO rreiated- 
Therefore, on large time scales P(9) can be approximated by a Gaussian which gives us first 
order term of R3. To calculate higher order terms which are not negligible on very long time 
scales the correlations between actions must also be taken into account. 

To further simplify the calculation of -R3, Ak{— rifce -7 ™' 1 ), in eq.(3.4), can be replaced 
by AiAjin^ 1 + rij 1 ) for the terms which survive due to n,/, = n\ + rij. Similarly for terms 
with rii = 7ij + rik or (rij = rij + n^), can be replaced by AiAj n 3 (l — rijn^ 1 ) and 
AjAj~ 1 rii(l — riinj 1 ) respectively. This leads us to following form of R%, 



here the 2 term corresponds to = rij + or = nj — rij, the 3 term corresponds to 
rij = rii + rik or rik = rij — rii . 

The Fourier transform of R% gives us the 3 rd order form factor K%, 



^3 = ^E^ 





(3.7) 
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K 3 (n, T 2 ) = J e 2 -[(^^3)r 1+ (r 2 -r3)r 2]jR3 ^ ^ ^ (g.g) 

= J e ^ir 1+ e 2 r 2 ] R3 ^ ^ d£idi2 ( 3 9 ) 

(where r\ — r 3 = £ x and r 2 — r 3 = £ 2 ). The eq.(3.9) follows from eq.(3.8) as R 3 depends only 
on differences r\ — r 3 , r 2 — r 3 . 

The further calculation of K 3 can be done by substituting eq.(3.7) in eq.(3.9) and by 
making use of following equalities (Appendix A) which follow from principle of uniformity, 

and 

i m v NJ g K J 

The result obtained depends on whether r's are greater or less than zero. This gives rise to 
following three possibilities. Case (1): both Ti,r 2 > or Ti,t 2 < 



-£(M + |r 2 |)e- 3T (3.13) 

It is obvious from above equation that K 3 falls very rapidly to zero for large T-values, that 
is, for time-scales on which sufficiently long periodic orbits exist in the phase space. On long 
time-scales, therefore this is similar to the RMT result (eq.(2.13)). Case (2): T\ > 0,r 2 < 

or T\ < 0, t 2 > 

( A h + 4% - - * (in - =f) « (n - =f ) (3.i4) 

^ 5 [|ri|/i(|r 2 |) + IralAdnl) - / 2 (|n|) - / 2 (|r 2 |)] e~ 3T (3.15) 

On substituting values of f\ = /(0, 1) and / 2 = /(0, 2) (Appendix A) in above equation, we 
get K 3 ~ (7 ( | ti | — |r 2 |)(e Ar72 — e JVri )e _3T where iV|ri| and N\t 2 \ are of the sam order as that 
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of T. This results into a nearly zero K 3 on large-time scales which is again similar to RMT 
result. 

Note the above-mentioned similarity between K 3 results for quantum maps and RMT 
has been shown here only for those time-scales at which principle of uniformity is well- 
applicable to the distribution of periodic orbits. No conclusion can be drawn about the 
short-time scales from the above analysis although the deviation of two-point fluctuation 
measures for quantum maps from those of RMT [2,7] suggests us to expect the same for 
higher orders too. 

B. 4 th -Order Correlation 

To calculate the 4 th order correlation function, We substitute eq.(2.4) and eq.(2.21) in 
eq.(2.3), with k = 4. This gives us 

i? 4 (Ei, E 2 , E 3 , E A ) = — ]T AiAjA k A r (exp [(m^ + rrijUj + m k n k + m r n r ) iE\) 

ijk mi,mj ,nik=±l 



exp 
exp 



(miUiii + mpifii + m k n k £ 3 ) — 



i' 

yTriiWi + rrijWi + m k W k + m r W r ) — 

h 



(3.16) 

Again the significant contributions to i? 4 come from following four types of terms, Case 



(1): Terms where pairwise cancellation occurs i.e terms with = n k ,nj = n r and Wi 
W k , Wj = W r (and their permutations) 



The contribution i? 4i from such terms can be written as follows, 

perm ij m=±l 

here J2 pe rm refers to the sum over all possible permutations of pairs. 



«« = ^EE A i A J E -P Wi - l 3 ) + n 3 l 2 )) (3.17) 



Case (2): Terms with rii — rij — n k — n r = (and other such permutations) 

The contributions to i? 4 from terms with — rij — n k — n r = can be written as follows, 
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a2 ( n2 i n k n j n l\ (2irmi 
E Af n 3 n k - -± 3 —± E ex P -^F~ n ^ - - n k £ 3 ) 



cxp ( ^ (Wi - ^ - W fc - W r ) ) ) (3.18) 



Similarly one can write the contributions from terms with rij — rii — n k — n r = and 
n k — ni — rij — n r = 0). The symbol R 4ii will refer to the sum of contibutions of all such 
terms. 

The contribution R 4 m from term with n r — rii — rij — n k — is 

■R'Ain = £ E ^A]A\ (— + - -) E (rtfi + n/ 2 + n k i 3 )) 

N \ n i n j n i n k n i n k J n ^ zl \ N J 

^exp (-™ (W r - Wi - W,- - W fc ))) (3.19) 
Case Terms with rii + rij — n k — n r = (and other such permutations) 



9 

N 4 



The contribution to R 4 from terms with n» + rij — n k — n r = can be written as follows, 
E A\A)n k (L + L-J*) E exp ( ^ + n/ 2 - n fe £ 3 )) 

(Wi + Wj-Wfc-Wr))) (3.20) 



' im 

Similarly one can write the contributions from terms with from terms rii — — n k + n r = 
and rii ~ n j + n k — n r = 0. The symbol R 4iv refers to the sum of contibutions of all such 
terms. 

Thus R 4 can be written as follows 

R 4 — Rii + Riii + Riiii + R&iv (3-21) 

where in each of the contributions R 4i ...R 4iv the terms of type < e v ^ +Wj+w/fc ~ w/r > can be 
replaced by e~ 4T (as done earlier for K 3 ). The Fourier transform of R 4 gives us the 4 rd order 
form factor K 4 , 

K 4 (ri, r 2 , r 3 ) = (2tt) 3 / e^^^^ ^ ^ d £ ld £ 2 d£ 3 ( 3 . 22 ) 
= + + ^ + tf iw (3.23) 
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where K, .„.,„.„■ = f e 2 ™^+^+ e ^R 4{miv) (4, 4a, 40 dhdi 2 d£ 3 . 

Now K 4 can be calculated by substituting eqs. (3. 17-3.20) in eq.(3.23) and using equalities 
eqs.(3.10,3.11). Again as for K 3 , the result depends on whether r's are greater or less than 
zero. This gives rise to following three possibilities. 

Case (1) Ti, r 2 , r 3 > or r u r 2 , r 3 < 

In this case, except for K 4iii , the contributions from all others, namely, K 4i , K 4ii and 
K 4iv are zero. Thus 




(3.24) 



^(M + N + |r 3 |K 4T (3.25) 

Due to presence of the exponentially decaying factor, K 4 turns out to be approximately zero 
for large T-values which is again similar to RMT results. 

Case (2) Any two oJti,t 2 ,t 3 positive (negative) and the third one negative (positive) 

Let Ti, Tj > and < (where k can take any of the values 1, 2, 3). In this case, 



K 4i ~g[5 (Ti + T k ) \Ti\\Tj\ + 5 (Tj + T k ) \Ti\\Tj\] (3.26) 

K<* = g HMrtihW) - A(M)/ 2 (N) - MN^dr,!)] e~ 4T (3.27) 

K 4iii = (3.28) 

K 4lv = g [NA(|r fe |) + \T \h(\T k \) - f 2 (\ Tk \)\ e~ 4T (3.29) 



As obvious from above equations, only K 4i does not contain an exponentially decaying factor 
and therefore makes a non-zero contribution to K 4 in large time-limits. The above results 
are valid also if n,Tj < and T k > 0. A comparison with RMT results shows that this 
lowest order contribution to K 4 is same as that in RMT for |ti| , |r 2 | , |r 3 | <^ 1. 
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C. /c th -Order Correlation 



The method used in calculation of 3 rd and A th order form factors can further be general- 
ized to the k th order correlation function The substitution of eqs.(2.4) and (2.21) in eq.(2.3) 
gives us 



n,..,ik 



mi ,..,mfc=±l 



(^exp 



2tt 



imi ( n h E + -j^n h ti + 



h 



....exp 



It can further be rearranged as follows, 

-,k k 



Rk 



4 e n-", e (-p 

il,..,ifci=l mi,..,m t =±l \ 

— ^min^i 



N k 
exp 



vm k ( n ifc £ + — 



(3.30) 



AT 



i=i 



[exp 



11 1=1 



(3.31) 



As due to ±1 values taken by each m h I = 1, 2, .., k there can be 2 k different combinations 
of ra i( 's in the 1 st exponent of eq.(3.31). Let M(r) be the set of a particular choice of values 
for each m ; in the set {mi,m 2 , ..m k } Therefore there can exist 2 k such sets, denoted by M(r) 
with r = 1 — > 2 fe , of which only 2 fc_1 sets are distinct. Here two sets M(r) and M(r') are 
considered indistinct if the values of each mi in M(r) is oppsite (in sign) to that in M(r'). 
Now eq.(3.31) can be rewritten as 

,fe 2 fc 



iv r=l 



(3.32) 



where Fjw r is given as follows, 



f~m t = Y.W A n ( ex p 

i\,..,i k l=l \ 
k-1 



iE ^ min k 
i=i 



exp 



.2tt 



J2 m l n i 



1=1 



exp 



ft 



(3.33) 



here, in above equation, the values taken by m/'s are the same as in M(r). 

Now for the same reason, as given for 3 rd order correlation function, only those terms 
of eq.(3.33) will contribute significantly to R k for which the multiplying factor of E in the 
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exponent is zero. Due to /c-summations over periodic orbits, each summation containing 
a large number of them, there are many possibilities, resulting in a zero coefficient of E. 
These various possibilities may arise, due to the "group- wise cancellation of periods" in 
the 1 st exponent of eq.(3.32), containing various groups of periods in the exponent, where in 
each group the positive traversals of few orbits are cancelled by the negative traversals of few 
other orbits. Let G stand for any division of indices 1, 2, .., k into q subgroups (Gi, G2, ••, G q ) 
then a term appearing in eq.(3.33) will make a non-zero contribution to Fm t if it satisfies 
following condition, 

J2mm H =0 (j = l,2,..g) (3.34) 

where the summation is over indices / present in subgroup Gj and so on. As obvious, one 
of these subgroups will contain index i^. Later we will need to distinguish it form other 
subgroups; let us call it G'y Thus F Mr can be rewritten as 

h 



7 M r = J2{ ex p 

G \ 



II C aj (3.35) 



1=1 

here J2g implies the summation over all possible divisions of indices ii,i2, --,ik into various 
subgroups, and flj implies the product of the contributions from all G/s for one such 
division, where 



h,-,3k \l=\ ' 1=1 



2171 V- K 



(3.36) 



here the p indices contained in subgroup Gj are denoted by ji,j2,--,jp with mi as the 
coefficient of rij l and bi's refer to £'s associated with these indices (e.g if ji = i\ then bi = t\ 
and so on); if ji = i k then b[ = bk = 0. Note that, in the group Gj, the index j 1 is 
always chosen to be ik (so as to simplify the presentation) with b\ = bk = 0. Now by using 



n 



— Z)f=2 a i n jn w hh ai — — m i/ m i, and Aj t ~ n^e 7 ™ J ', one can show that 

1=1 1=2 u ji s=2^l IL is 

By using eq.(3.37), Cgj can further be reduced as follows 

20 



' A ai + 1 p A a s+ 1 

e ie«.-&t n h 



n 



s=2^l H js 



exp 



2m p 

^r(E( 6 ' -bi)m injl ) 

iV 1=2 



(3.38) 



h,-,ji \l=2 -"ji 

Note here that for j — f, b\ = 0. It can be seen from above equations that the most 
significant contribution to R k comes from those terms where a pairwise cancellation of 
time-periods rij, rij, rik, n s , ... (appearing as a factor of E in the exponent) as well as actions 
Wi, Wj, Wk, W s , ... (so that there is no exponential decay) occurs which is possible only when 
k is even. The contribution R ki from such terms can be written as follows, 

,k 



9 

N k 



E E A l A l : 

P erm ji,-,jk 

\k-2)/2 

E n m-i ( b 2i-i ~ 6a) + Wjfc-i&k-i 



E ex P 

m=±l 



2nmi 



N 



i=i 



(3.39) 



The I] perm implies the summation over all possible permutations of indices j±, j 2 , ..,jk taken 
from set ii,i 2 ,--,ik- The contributions from all other terms contain an exponential term, 
with sum over actions (in units of H), as its exponent (i.e term < exp(X)raj.W^) >• In 
large time limits, the application of CLT again permits us to replace this term by e~ kT . The 
contribution of such terms to Rk can be given as follows 

Mr G 3 ( P 



k Mr u j I P A ai+l p A 01 ^ 1 

^=4eeii e E««zS=r n 1 



exp 



r G j ji,..,j p \l=2 

'2m t p 



n 



a l - 
31 



N 



\l=2 



s=2^l U js 



-kT 



(3.40) 



The substitution of eqs.(3.39) and (3.40) in eq.(2.9) gives us following result for K\. 



K kl + K, 



kii 



(3.41) 



where 



and 



K 



hi 



N k 



E E A h A h— A 3k-i 



P erm jl,-,jk 



E <K*2Jfc-l 
m=±l 



fc-i 



n - 

i=i 



N 



)6(t 2 l-l + T 2 l) 



(3.42) 
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' r G j h,..,j p \l=2 "j. s=2^l a js 



M r 



Gj 



v 

w 

U=2 



n + 



AT ■ 



-fer 



(3.43) 



where t/ is the r- variable associated with 6/ with ti = if bi = t k . By using eq.(A3), one can 
further reduce above equations in following form, 



K ki = J2i 6 ( t pi + T P2) S ( r P 3 + t pJ -(k/2 - l)terms] | 

T Pi I \ T P3 I ••• \ T P(k-l) \$k, 



(3.44) 



where J2p refers to sum over all possible permutations of indices {pi,P2, --,^4} taken from 
set {1, 2, .., k — 1}, and 



M(r) 

K kii = gY^ Y[D Gj e 

r G j 



-kT 



(3.45) 



with Dq. given as follows, 



H 6 (hi + mi) 



1=1 U=2 
' P 

Y[S(hi - mi) 

.1=2 



5>,/(a, + l,a,-l) 1] /(«. + !,«.) ^ ^ ^ (3-46) 



1=2 



s=2,# 



J2aif(ai + l,ai - 1) J] /(a, + 1, a s ) Gj = Gj 



(3.47) 



Z=2 s=2,^l 

with /i; = sign(ti). Here, as obvious, the contribution K ki exist only for for all even order 
form factors. For odd order form factors, only K kii contributes. But, on large time scales, the 
presence of an exponential decaying factor makes this contribution very small. Thus, on large 
time-scales (i.e n — > 00, iV — > 00, n/N < 1) one can write 



K. 



k—odd 



~ 



(3.48) 



A comparison of eqs.(3.44) and (3.47) with eqs.(2.14) and (2.13), informs us that the results 
obtained for both odd as well as even order form factors agree well with that of RMT. 



IV. NUMERICAL VERIFICATION 

This section contains a numerical study of 3 rd and 4 th order fluctuation measures, namely, 
skewness 71 and excess 72. 
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We choose the kicked rotor system for this purpose as it has been an active model of 
research, containing a variety of features such as localization, resonance, dependence of the 
spectra on number theoretical properties etc and has been used as a model for a very wide 
range of physical systems. For a better understanding, we briefly review the quantum and 
classical mechanics of kicked rotor in this section. 



A. The Kicked Rotor: Classical and Quantum Dynamics 



The kicked rotor can be described as a pendulum subjected to periodic kicks (with period 
T) with the following Hamiltonian 

(p + if 



H 



+ K cos(6 + e ) J2 S(t-nT) 



(4.1) 



where K is the stochasticity parameter. The parameter 7 and 9 are introduced in the 
Hamiltonian in order to mimick the effects of the time-reversal (T) and the parity (P) 
symmetry breaking in the quantum system. 

The related quantum dynamics can be described, by using Floquet's theorem, by a 
discrete time evolution operator U = BG where B = exp(—iKcos(9 + 9 )/k) and G = 
exp(—i(p + 7) 2 /4fr). The nature of the quantum dynamics and therefore the statistical 
properties of the associated quantum operators depend on h and K. For a rational value of 
hT/2n, the dynamics can be confined to a torus while for irrational value it takes place on 
a cylinder. We employ torus boundary conditions (q' — q + 2n,p' = p + 2nM/T) by taking 
hT/2Tc = M/N; both p and 9 then have discrete eigenvalues and U can be reduced to a 
finite iV-dimensional matrix of the form [32-34] 



U mn — ^.exp 



K (2i\m n 

. tT cos[— + eo 



exp 



l(m — n) 



N 



(4.2) 



(n 2 M 2 - vr 7 /)" 

where Ni — (N — l)/2 (with N odd) and m,n = -N 1: -JVi + 1, iV x . 

The quantum dynamics has a time reversal symmetry T for 7 = and a parity symmetry 
P for #0 = 0. Though the T-symmetry may be violated for 7 7^ 0, still a more generalized 
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antiunitary symmetry S = TP = PT can be preserved in the system if 9 — [14,2]. For 
K 2 » Nh, the quantum dynamics is delocalized in the phase space and two point spectral 
fluctuation measures have been shown to be well modeled by various symmetry classes of 
RMT [14,2]. In the opposite limit of weak chaos, namely, K 2 « Nh, the eigenstates 
localize in the momentum space and one obtains a Poisson distribution for the spectrum [2] . 
For numerical comparison of higher order measures with RMT, therefore, we choose various 
parameters such that condition K 2 » Nh is always satisfied. 



B. Numerical Study of Skewness and Excess 

In this section, we numerically study skewness and excess, that is, the 3 rd and 4 th order 
spectral fluctuation measure, for QKR-spectra and compare them with corresponding RMT 
results. Although the k th order form factor can easily be calculated for quantum maps (see 
eq.(2.10)), we choose to study 7! and 72 as corresponding numerical results (for form factors) 
in RMT are not available. But as both 7x and 7 2 can analytically be expressed in terms of 
3 rd and 4^-order correlations _R 3 and _R 4 (see ref.[15] for these expressions) and, therefore, 
are related to form factors K 3 and K 4 (eq.2.8), any conclusion about the validity of RM 
model for the former measure will be applicable for the latter too. 

Both skewness and excess are the functions of 3 rd and 4 th order central moment of a 
distribution respectively which contain the information about the probability of higher order 
events as compared to the lower order. More precisely, the skewness denotes the absence of 
symmetry in the distribution and can be defined as follows, 

TlM = ^ (4-3) 

The excess 72 describes the difference between the kurtosis values (i.e the fourth central 
moment calculated in units of the square of second central moment) of the distribution and 
that of a normal distribution (kurtosis=3), 
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where a 2 is the variance in the number of levels in a length of r mean spacings and /X3 and 
/i4 are corresponding 3 rd and 4 th central moments. If excess is less than zero, the curve is 
platykurtic and if it is positive, leptokurtic. 

For both 71 and 72 studies, the spectral data consists of the eigenvalues of 50 matrices 
of dimension N = 199 obtained by diagonalizing [/-matrix (eq.(4.2)) for various values of 
K in the neighbourhood of K ~ 20000. The choice of such a high value of K is made to 
ensure the derealization of quantum dynamics which makes U a full random matrix (see 
ref.[2]). Due to strong sensitivity of the eigenvalues to small changes in K, these sequences 
of quasienergies can be regarded as mutually independent. 

The results obtained for 7! and 72 of the QKR-spectra are displayed, as function of r 
in figures 1 and 2;. For comparison, the corresponding RMT-results (taken from ref. [15]) 
are also given in each of the figures. The good agreement between QKR and RMT results, 
indicated by each of these figures reconfirms our semiclassical results obtained in sectionpl. 
Furthermore, as indicated by these figures, the higher order correlations e.g 3 rd and 4 th 
order, seem to be very weak, nearly zero, at very long energy ranges and levels seem to be 
uncorrelated while strong correlation seem to be existing for short ranges. A continuously 
decreasing value of 71 implies the tendency of the distribution to appear more and more 
symmetric as range of the distribution increases and finally acquire the Gaussian form for 
long energy ranges. This conclusion is also supported by our 72 study which shows a larger 
probability of higher order events for QKR at small energy ranges (as compared to Gaussian 
case) while, for large energy ranges, it acquires the same form for both. 

V. CONCLUSION 

We conclude this article with a summary of our principle results and a brief discussion 
of the open problems. 

We have shown that, in long time limits, the higher order spectral correlations in quan- 
tum chaotic maps can be well modeled by corresponding ones in RMT. By using the example 
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of kicked rotor, we have verified this numerically too. These results are also valid for conser- 
vative Hamiltonians; this follows due to the similar expressions for the level density in terms 
of classical periodic orbits for both the cases (see Ref.[4] for level density for autonomous 
case). The various summation formulas used to evaluate fluctuation measures such as sum 
over amplitudes (Appendix A), whose derivation depends on the uniform distribution of 
periodic orbits and therefore only on the strongly chaotic nature of dynamics, still remain 
valid. Thus one obtains the similar results for form factors in autonomous systems. This 
should not be surprising, as both Gaussian ensembles (the ensembles of hermitian matrices 
and therefore of conservative hamiltonians) and circular ensembles are known to have the 
same statistical behaviour in semiclassical limit [16]. 

Although the study presented here deals with higher order correlations at fixed value of 
parameter, we expect the validity of random matrix model also for higher order parametric 
correlations in quantum chaotic spectra, but again only on long time-scales. The intuition, 
based on the analogy of 2 nd order parametric density correlations for mesoscopic systems 
with a disordered potential, one dimensional many body hamiltonian (Calogero-Sutherland 
model), scattering systems, quantum chaotic systems and random matrix models, as well 
as the existence of a common mathematical base (nonlinear a-model and supersymmetry 
approach) further encourages us to hope the extension of this analogy for higher order 
correlations too. 

The ignorance about action correlations on very long time scales handicaps us from doing 
the same analysis on these scales. It will also be of interest to study these correlations on 
short time scales. It is on these scales where 2 nd order spectral correlations deviate from 
RMT and show non-universal behaviour. One expects to see similar deviations for higher 
orders too. For a complete understanding of higer order correlations in quantum chaotic 
systems, therefore, one should study the action and periodic orbit correlations; this will give 
us some insight in the behaviour of fluctuations in other systems (above-mentioned) too. 
We intend to do so in future. 
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APPENDIX A: CALCULATION OF F(A, B) 
To evaluate the following sum 

/M) = e3'(M-£) (AD 

where A, is the amplitude and rij is the period of the orbit in the classical phase space. We 
remind ourselves that, in the phase space, the orbits proliferate exponentially with time and 
the density of distribution of periods over long orbits is given by exp(27|n|)/|n| with 7 as 
the entropy of the classical motion. In large time-limits. This, along with the approximation 
Aj ~ rijexp{— 7 1 7^- 1), enables us to make the following replacement 

A a n ■ 00 n 

£ ^ (|T| " N } ~ £ e M-m(a - 2))n a - b - 1 5(\r\ - -) (A2) 

j n j JV n=0 iV 

which gives 

f(a, b) = exp(- 7 iV|r|(a - l))\r\ M N^ b (A3) 
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FIGURES 

FIG. 1. The behaviour of 71 (r) with respect to r, with N = 199, h = 1, T = 1, 
X = 20002 -» 20050 and for (a) 7 = 0.0,00 = tt/2N, TR-preserved (b) 7 = 0.7071, = tt/2N, 
TR-broken The solid curve depicts the corresponding RMT behaviour, namely COE (in 1(a)) and 
CUE limit (in 1(b)). 

FIG. 2. The behaviour of 72 (r) with respect to r, with N = 199, h = 1, T = 1, 
K = 20002 -» 20050 and for (a) 7 = 0.0,6*0 = tt/ 2 ^ TR-preserved (b) 7 = 0.7071, O = vr/2iV, 
TR-broken The solid curve depict the corresponding RMT behaviour, namely COE (in 2(a)) and 
CUE limit (in 2(b)). 
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